Impacts of Groundwater Pumping for Hydraulic Fracturing on Aquifers Overlying the Eagle Ford Shale.

Hydraulic fracturing (HF) events consume high volumes of water over a short time. When groundwater is the source, the additional pumping by rig/frack supply wells (RFSWs) may impose costs on owners of other sector wells (OSWs) by lowering the hydraulic head. The Carrizo-Wilcox aquifer in south Texas is the main source of water for HF of the Eagle Ford Shale (EFS) Play. The objectives are to assess the impacts of groundwater pumping for HF supply on: (1) hydraulic heads in OSWs located nearby an RFSW and (2) volumetric fluxes between layers of the regional aquifer system compared to a baseline model without the effect of RFSW pumping. The study area spans the footprint of the EFS Play in Texas and extends from 2011 to 2020. The pumping schedules of 2500 RFSWs were estimated from reported pumped water volumes to supply 22,500 HF events. Median annual drawdowns in OSWs ranged from 0.2 to 6.6 m, whereas 95th percentile annual drawdowns exceeded 20 m. The magnitudes of drawdown increased from 2011 to 2020. Of the four layers that comprise the Carrizo-Wilcox aquifer, the upper Wilcox was the most intensively pumped for HF supply. During the peak HF year of 2014, the net flux to the upper Wilcox was 292 Mm3 compared to the baseline net flux for the same year of 278 Mm3 -a relative gain of 14 Mm3 . Pumping for HF supply has the potential to negatively impact nearby OSWs by capturing water from adjacent aquifer layers.

Optimizing Multiwell Aquifer Storage and Recovery Systems for Energy Use and Recovery Efficiency.

As more aquifer storage and recovery (ASR) systems are employed for management of water resources, the skillful operation of multiwell ASR systems has become very important to improve their performance. In this study, we developed MODFLOW and MT3DMS models to simulate a multiwell ASR system in a synthetic aquifer to assess effects of hydrogeological and operational factors on the performance of the multiwell ASR system. We evaluated a simplified (dual well) ASR system in comparison with complex system (three-, four-, five-, and seven-well systems). Recovery and energy efficiencies were calculated using the model simulations. Factors such as higher hydraulic conductivity and longitudinal dispersivity significantly reduced the recovery and energy efficiencies of the system. In contrast, increasing the volume of recharged water increased the recovery efficiency; however, the energy efficiency was reduced. Recovery and energy efficiencies also plummet when there is an increase in the underlying regional gradient and the designed storage duration. Operating the system multiple times can yield higher volume of potable water, but the energy efficiency may not vary significantly after the second operating cycle. Single-well systems and multiwell systems exhibit similar responses to changes in physical factors, although operational factors have a more pronounced effect on the multiwell systems. One of the major findings was that fewer wells in a multiwell ASR system can yield higher volume of potable water and better output with respect to the electrical power being consumed. The results provide design engineers with guidelines for optimizing performance of the multiwell ASR systems.

Precipitation mediates sap flux sensitivity to evaporative demand in the neotropics.

Transpiration in humid tropical forests modulates the global water cycle and is a key driver of climate regulation. Yet, our understanding of how tropical trees regulate sap flux in response to climate variability remains elusive. With a progressively warming climate, atmospheric evaporative demand [i.e., vapor pressure deficit (VPD)] will be increasingly important for plant functioning, becoming the major control of plant water use in the twenty-first century. Using measurements in 34 tree species at seven sites across a precipitation gradient in the neotropics, we determined how the maximum sap flux velocity (vmax) and the VPD threshold at which vmax is reached (VPDmax) vary with precipitation regime [mean annual precipitation (MAP); seasonal drought intensity (PDRY)] and two functional traits related to foliar and wood economics spectra [leaf mass per area (LMA); wood specific gravity (WSG)]. We show that, even though vmax is highly variable within sites, it follows a negative trend in response to increasing MAP and PDRY across sites. LMA and WSG exerted little effect on vmax and VPDmax, suggesting that these widely used functional traits provide limited explanatory power of dynamic plant responses to environmental variation within hyper-diverse forests. This study demonstrates that long-term precipitation plays an important role in the sap flux response of humid tropical forests to VPD. Our findings suggest that under higher evaporative demand, trees growing in wetter environments in humid tropical regions may be subjected to reduced water exchange with the atmosphere relative to trees growing in drier climates.

Leaf surface traits and water storage retention affect photosynthetic responses to leaf surface wetness among wet tropical forest and semiarid savanna plants.

While it is reasonable to predict that photosynthetic rates are inhibited while leaves are wet, leaf gas exchange measurements during wet conditions are challenging to obtain due to equipment limitations and the complexity of canopy-atmosphere interactions in forested environments. Thus, the objective of this study was to evaluate responses of seven tropical and three semiarid savanna plant species to simulated leaf wetness and test the hypotheses that (i) leaf wetness reduces photosynthetic rates (Anet), (ii) leaf traits explain different responses among species and (iii) leaves from wet environments are better adapted for wet leaf conditions than those from drier environments. The two sites were a tropical rainforest in northern Costa Rica with ~4200 mm annual rainfall and a savanna in central Texas with ~1100 mm. Gas exchange measurements were collected under dry and wet conditions on five sun-exposed leaf replicates from each species. Additional measurements included leaf wetness duration and stomatal density. We found that Anet responses varied greatly among species, but all plants maintained a baseline of activity under wet leaf conditions, suggesting that abaxial leaf Anet was a significant percentage of total leaf Anet for amphistomatous species. Among tropical species, Anet responses immediately after wetting ranged from -31% (Senna alata (L.) Roxb.) to +21% (Zamia skinneri Warsz. Ex. A. Dietr.), while all savanna species declined (up to -48%). After 10 min of drying, most species recovered Anet towards the observed status prior to wetting or surpassed it, with the exception of Quercus stellata Wangenh., a savanna species, which remained 13% below Anet dry. The combination of leaf wetness duration and leaf traits, such as stomatal density, trichomes or wax, most likely influenced Anet responses positively or negatively. There was also overlap between leaf traits and Anet responses of savanna and tropical plants. It is possible that these species converge on a relatively conservative response to wetness, each for divergent purposes (cooling, avoiding stomatal occlusion, or by several unique means of rapid drying). A better understanding of leaf wetness inhibiting photosynthesis is vital for accurate modeling of growth in forested environments; however, species adapted for wet environments may possess compensatory traits that mitigate these effects.

Mapping potential groundwater-dependent ecosystems for sustainable management.

Ecosystems which rely on either the surface expression or subsurface presence of groundwater are known as groundwater-dependent ecosystems (GDEs). A comprehensive inventory of GDE locations at an appropriate management scale is a necessary first-step for sustainable management of supporting aquifers; however, this information is unavailable for most areas of concern. To address this gap, this study created a two-step algorithm which analyzed existing geospatial and remote sensing data to identify potential GDEs at both state/province and aquifer/basin scales. At the state/province scale, a geospatial information system (GIS) database was constructed for Texas, including climate, topography, hydrology, and ecology data. From these data, a GDE index was calculated, which combined vegetative and hydrological indicators. The results indicated that central Texas, particularly the Edwards Aquifer region, had highest potential to host GDEs. Next, an aquifer/basin scale remote sensing-based algorithm was created to provide more detailed maps of GDEs in the Edwards Aquifer region. This algorithm used Landsat ETM+ and MODIS images to track the changes of NDVI for each vegetation pixel. The NDVI dynamics were used to identify the vegetation with high potential to use groundwater--such plants remain high NDVI during extended dry periods and also exhibit low seasonal and inter-annual NDVI changes between dry and wet seasons/years. The results indicated that 8% of natural vegetation was very likely using groundwater. Of the potential GDEs identified, 75% were located on shallow soil averaging 45 cm in depth. The dominant GDE species were live oak, ashe juniper, and mesquite.

A statistical method for estimating wood thermal diffusivity and probe geometry using in situ heat response curves from sap flow measurements.

The heat pulse method is widely used to measure water flux through plants; it works by using the speed at which a heat pulse is propagated through the system to infer the velocity of water through a porous medium. No systematic, non-destructive calibration procedure exists to determine the site-specific parameters necessary for calculating sap velocity, e.g., wood thermal diffusivity and probe spacing. Such parameter calibration is crucial to obtain the correct transpiration flux density from the sap flow measurements at the plant scale and subsequently to upscale tree-level water fluxes to canopy and landscape scales. The purpose of this study is to present a statistical framework for sampling and simultaneously estimating the tree's thermal diffusivity and probe spacing from in situ heat response curves collected by the implanted probes of a heat ratio measurement device. Conditioned on the time traces of wood temperature following a heat pulse, the parameters are inferred using a Bayesian inversion technique, based on the Markov chain Monte Carlo sampling method. The primary advantage of the proposed methodology is that it does not require knowledge of probe spacing or any further intrusive sampling of sapwood. The Bayesian framework also enables direct quantification of uncertainty in estimated sap flow velocity. Experiments using synthetic data show that repeated tests using the same apparatus are essential for obtaining reliable and accurate solutions. When applied to field conditions, these tests can be obtained in different seasons and can be automated using the existing data logging system. Empirical factors are introduced to account for the influence of non-ideal probe geometry on the estimation of heat pulse velocity, and are estimated in this study as well. The proposed methodology may be tested for its applicability to realistic field conditions, with an ultimate goal of calibrating heat ratio sap flow systems in practical applications.

Adaptive evolution of plastron shape in emydine turtles.

Morphology reflects ecological pressures, phylogeny, and genetic and biophysical constraints. Disentangling their influence is fundamental to understanding selection and trait evolution. Here, we assess the contributions of function, phylogeny, and habitat to patterns of plastron (ventral shell) shape variation in emydine turtles. We quantify shape variation using geometric morphometrics, and determine the influence of several variables on shape using path analysis. Factors influencing plastron shape variation are similar between emydine turtles and the more inclusive Testudinoidea. We evaluate the fit of various evolutionary models to the shape data to investigate the selective landscape responsible for the observed morphological patterns. The presence of a hinge on the plastron accounts for most morphological variance, but phylogeny and habitat also correlate with shape. The distribution of shape variance across emydine phylogeny is most consistent with an evolutionary model containing two adaptive zones--one for turtles with kinetic plastra, and one for turtles with rigid plastra. Models with more complex adaptive landscapes often fit the data only as well as the null model (purely stochastic evolution). The adaptive landscape of plastron shape in Emydinae may be relatively simple because plastral kinesis imposes overriding mechanical constraints on the evolution of form.

Modeling vadose zone processes during land application of food-processing waste water in California's Central Valley.

Land application of food-processing waste water occurs throughout California's Central Valley and may be degrading local ground water quality, primarily by increasing salinity and nitrogen levels. Natural attenuation is considered a treatment strategy for the waste, which often contains elevated levels of easily degradable organic carbon. Several key biogeochemical processes in the vadose zone alter the characteristics of the waste water before it reaches the ground water table, including microbial degradation, crop nutrient uptake, mineral precipitation, and ion exchange. This study used a process-based, multi-component reactive flow and transport model (MIN3P) to numerically simulate waste water migration in the vadose zone and to estimate its attenuation capacity. To address the high variability in site conditions and waste-stream characteristics, four food-processing industries were coupled with three site scenarios to simulate a range of land application outcomes. The simulations estimated that typically between 30 and 150% of the salt loading to the land surface reaches the ground water, resulting in dissolved solids concentrations up to sixteen times larger than the 500 mg L(-1) water quality objective. Site conditions, namely the ratio of hydraulic conductivity to the application rate, strongly influenced the amount of nitrate reaching the ground water, which ranged from zero to nine times the total loading applied. Rock-water interaction and nitrification explain salt and nitrate concentrations that exceed the levels present in the waste water. While source control remains the only method to prevent ground water degradation from saline wastes, proper site selection and waste application methods can reduce the risk of ground water degradation from nitrogen compounds.